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Abstract 

We present higher order polynomial algebras which are the dynamical symmetry 
algebras of a wide class of multi-mode boson systems in non-linear optics and laser 
physics. We construct their unitary representations and the corresponding single- 
variable differential operator realizations. We then use the results to obtain exact 
(Bethe ansatz) solutions to the multi-mode boson systems, including the Bose- 
Einstein Condensate models as the special cases. 
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1 Introduction 

This paper is the second of a series of three devoting to polynomial algebraic structures 
and exact solutions of general quantum non-linear optical models. In the previous work 
[1], we exactly solved a class of two-mode boson systems by identifying and applying the 
underlying dynamical polynomial algebra symmetries. In this article, we generalize the 
results of [1] and deal with a general class of multi-mode boson systems. 

There has been constant interest in the study of deformed Lie algebras and their 
applications to physics. This has been motivated somewhat by the realization that the 
traditional linear Lie algebras and their affine counterparts may be too restrictive as there 
is no reason for symmetries to be linear in the first place. 

Polynomial algebras are a particular type of non-linear deformations of Lie algebras [2] . 
Their applications span across a diverse field of theoretical physics whereby they appear 
in topics such as quantum mechanics, Yang-Mills type gauge theories, quantum non- 
linear optics, integrable systems and (quasi-)exactly solvable models, to name a few (see 
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e.g. [3]-[13]). Different realizations of polynomial algebras have been extensively studied 
in [10, 11, 12] (see also [14, 15, 16, 17] for differential operator realizations of certain 
quadratic and cubic algebras in connection with quasi-exact integrability [18, 19, 20]). 

One of the main results of this paper is the derivation of the exact eigenfunctions and 
energy eigenvalues of the multi-boson Hamiltonian, 

r+s r+s 
i i<j 

where and throughout Oj (a|) and — a\ai are bosonic annihilation (creation) and num- 
ber operators, respectively, and Wi, Wij and g are real coupling constants. Hamiltonians 
of the form (1.1) appear in the description of various physical systems of interest such 
as non-linear optics and laser physics. Using the Bargmann representation, it was shown 
in [21] that the multi-boson system is quasi-exactly solvable (see also [22]). Some special 
cases of the system have been studied in e.g. [23, 24] by means of the Algebraic Bethe 
Ansatz method. Here we determine the underlying dynamical polynomial algebra symme- 
tries and the exact solutions of the general Hamiltonian (1.1). The polynomial algebras 
and representations constructed in this paper are of interest in their own right, in view of 
their potential applications in other fields. 

This paper is organized as follows. In sections 2 we introduce higher order polynomial 
algebras with multi-mode boson realizations. We construct their unitary representations 
in the Fock space and the corresponding single- variable differential operator realizations. 
We then identify the polynomial algebras as the dynamical symmetries of Hamiltonian 
(1.1) in section 3, and solve for the eigenvalue problem via the Functional Bethe Ansatz 
method (see e.g. [25, 26, 27]). In section 4, we present explicit results for models cor- 
responding to the special cases r -I- s < 4, providing a unified treatment of the so-called 
Bose-Einstein Condensate (BEC) models. We summarize our results in section 5 and 
discuss further avenues of investigation. 



2 Polynomial deformations of sl(2) algebra 

We first briefiy review the polynomial algebras defined in [1] . Let A; be a positive integer, 
A; = 1, 2, • • •. Consider the polynomial algebra of degree k — 1 defined by 

[Qo,Q±] = ±Q±, 

[Q+,Q-] = <P^'\Qo)-<P^'\Qo-i), (2.1) 

where 

.^"=' «2o) = - n (oo + ^) + n - ^) (2.2) 



2 



is a polynomial in Qq of degree k. The algebra admits Casimir operator of the form 

C = Q.Q+ + 0(')(go) = Q+Q- + - 1). (2.3) 

As was shown in [1], algebra (2.1) has an infinite dimensional irreducible unitary repre- 
sentation given by the following one-mode boson realization, 

«-^(^'"**'' "^-^w*"''' p-^' 

In this realization, the Casimir (2.3) takes fixed value 

c=n(^-^). (2.5) 



which corresponds to the common eigenvalue [1] for the direct sum of k irreducible repre- 
sentations with quantum number q — • • • , (*^~^)*^+^ in each case. Note that for 
k — 2, C becomes 3/16 and q reduces to the well-known quantum number K — 1/4, 3/4 
used in su{l, 1). Therefore quantum number g is a very natural choice for labeling the 
Fock states, denoted as |g, n), of the irreducible representations of (2.1). Explicitly, 

hn)^ 10), n = 0,l,---. (2.6) 

^/[k{n + q-^)]\ 

The action of Qq, Q± on these states are as follows 

Qo\q,n) = iq + n)\q,n), 

^ , V A / ik-lV^ . 



i=l 



Q-M = l[(n + q-^^-^^t^y \q,n-l). (2.7) 



i=l 



Now let us take r mutually commuting copies of algebra (2.1), {Q+ ,Q-\Qo^}, of 
degree ki — 1, with the multi-boson realization {i — 1, ■ ■ ■ ,r) 

Consider the algebra generated by 

Q+ = gi^)---gi:\ Q_ = g«---gL^\ Qo = -EQ?. (2.9) 



r i=i 



It can be shown that Qo,± form the polynomial algebra of degree {J2i=i h ~ 

[Qo,Q±] = ±Q±, 

= -n0(^^)(Qo,{£}) + n0('=^)(Qo-l,{>C}). (2.10) 
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In the above, {£} stands for the set {£i, • • • , Cr-i}, with 

A- = Q^^-Qt'^ (2.11) 

being the r — 1 central elements of algebra (2.10), [Cj, Qo,±] = 0, and 

<P^'HQo, {c}) = fl f 2o + E - 1 E /X/:. + (2.12) 

is a polynomial in Qq and >Cj of degree ki. Prom (2.11), we can show that 

= E A- + Qi^^ , Qi^^ = Qo - - E ^'^^r (2-13) 

j=i ^ j=l 

Unitary irreducible representations of polynomial algebra (2.10) are infinite dimen- 
sional. The corresponding Fock states take the form HLi where rii — 0, 1, • • • , 
and Qi — • • • , feii^ili. By means of the central elements (2.11), we can show 
that the Fock states become 

^ki(n+qr+Si-^) 

\{Qh {I}) = n s, + qr-qi + n)^f[ , ' , , (2.14) 

i=l i=l Y [Ki{n -\- Qr -t Si — p-JJ ! 

where n = 0, 1, ■ ■ ■, and Si = J^jZl Ij for i = 1,2, ■■■ ,r ~ 1 and Sr = 0; {q} = {qi, - ■ ■ , qr} 
and {/} = {^1, • . . , Ir-i} are the constant values taken by the central elements £j. Note 
that Si + qr — qi > (so qr + Si - > 0). 

i 

We now construct polynomial algebra which has finite dimensional unitary irreducible 
representations and thus can be identified as the dynamical symmetry algebra of the multi- 
boson Hamiltonian (1.1). We do so by considering two mutually commuting algebras 
{Qi%} of degree ELi h - 1 and {Q^'i} of degree ^llr+i ' 1' ^i^h 

0(i)_-rT (4)^' oW-TtJ^ oW-ix-l/'iV. + l 

( r,'^\ki r+s ( ^ \ki 1 r+s i / -i 

i=r+l iVh)''^ i=r+l (V^)*' « ii^^^l ^ 

Introducing new generators. 



(2.15) 



P+ = QS'qL^ P_ = qJ)Q«, Po = ^K^-QfO- (2.16) 

We can easily show that Vq^± form a polynomial algebra of degree (X^i^i^^^ — 1) which 
close under the following commutation relations: 

[p+^^-\ = -riv^S''\^o,/c,{£}W) n%r)(Po,/c,{£}(2)) 

i=\ i=r+\ 

+ i{^f\Vo-l,lC,{CY^^) n ^t\Vo-l,K.,{CY^^), (2.17) 

i=l i=r+l 
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where {C}^^^ = {£i, • • • , JCr-i} and {Cy^^ = {Cr+i, • • • , jCr+s-i} with >Cj given by (2.11) 
are the central elements of (2.17), i.e. [Cj,Vo^±] — 0, 

/C = ^ + (2.18) 
is the (r + s — l)-the central element of the algebra, [/C, Po,±] = 0, and finally 

ki / r-1 1 r-1 -7. 















ill 





^ r ^ ' ^ ' 1.2 

lj,=i ' M=l « 

r+s— 1 1 r+s— 1 



«i / r+s-i 1 r+s-i -7, 



(2.19) 



Unitary irreducible representations of polynomial algebra (2.17) arc finite dimensional. 
To show this, let \{qY'^\n^^\ {/}(^)) and \{q}^^\n^^\ {1^^^) be the Fock states of the alge- 
bras |Qq^j_| and |Qo^±| respectively, where nS^\ n*-^^ = 0, 1, ■ ■ ■, and {g}*-^-* = {qi, ■ ■ ■ , qr}, 
= {qr+i, ■ ■ ■ , qr+s}- The representations of {Vo^±} arc then given by the Fock states 
\{q}'^^\ n'^^\ {l}^^^) \{q}'''^\ n'^'^\ {O*"^"*)- Since /C is a central element of the algebra, it must 
be a constant, denoted as k below, on any irreducible representations. This imposes the 
following constraint 

^(1) + ^(2) ^2K-qr- qr+s - t, (2.20) 



where 



-1 r -| r+s 



.(1) 



^ i=l * i=r+l 

^ = l,---,r-l, s« = 0, 

j=i 



4'^ = ^.-^ ^ = r + l,---,r + s-l, 41 = 0. (2.21) 

Let J\f = 2k, — qr — qr+s — t. Then obviously A/" = 0, 1, 2, • • • , and the Fock states become 



\{qY'\{q}^'\n,{lY'\{lY'\K) = n, \{qY'\ J\f - n, {lY'^) 

nLi(al) ' ^ Y\^+M) ^ ' ^ >) 

^nLi[^i(n + + - ^)]! nii%i[^i(2/^ -qr-t + sf^ - ^ - n)]! ' 

(2.22) 

where n = 0, 1, • • • , jV. This gives us the Af +1 dimensional irreducible representation of 
(2.17), 

vm^'\{<iY'\n.{iY'\{iY'\^) 
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V^\{qY'\{qY'\n,{lY'\{lY'\K) 



j=r+li=l \ '^j / 



X n n h + + + i^^>^'^' -f^}^'^' ^ + 1' -f^}^'^' -f^}^'^' '^)' 

,7 = 1 j=l V % / 



(2) , ^l^j 



iki — 1 



r+s kj / 

n \{[2n-n-qr-t-rsf-r 



X n n \n + qr + s'P - p n - 1, {/}«, {/}(^), «:). 

(2.23) 



By using the standard Fock-Bargmann correspondence, we can map the Fock states 
{qY^\n, {lY^\ k) to the monomials mz = Ul=i z^' / 111=^1 ^i'^ 

^JUUMn + qr + s\'^ - ^)]! nr+Vi[^i(2« - - i + - -| - n)]! 

n = 0, 1, • • • , A/". (2.24) 
The corresponding single- variable differential operator reahzation of (2.17) is 

d 1 (1) 

Vo = Z— - K + qr + -2_^S)', 

az J, ^ J 



-^5=^ n li f^::^ + + 4' - '^^) (-5) 



ns^, (v^ A *J 



These differential operators form the same (A/" + 1) dimensional representations in the 
space of monomials as those reahzed by (2.15) in the corresponding Fock space. We 

remark that because 11^=1 IliLi (^Qr + ■s^-^'' — = ^r all the allowed values 

(noting s^^^ = 0) there is no z^'^ term in above and thus the differential operator 
expressions (2.25) are non-singular. 
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3 Exact solution of the multi-mode boson systems 

We now use the differential operator realization (2.25) to exactly solve the multi-mode 
boson Hamiltonian (1.1). 

By (2.16) and the multi-boson realization (2.15), we may express the Hamiltonian 
(1.1) in terms of the generators of the polynomial algebra (2.17), 

H = Y,w^N, + E^^u^^^i + ^ n (yS)'^ (^+ + ^-) (3-1) 

i i<j i=l ^ ' 

with the number operators having the following expressions in Pq and L 

^r— 1 p 1 v^r— 1 ,, r \ 1 



^ \ A:, /C + n + Epi/:.-^Epl/x^.j-i for. = l,...,r 
' I A;,(/C-Po + E;±r'>C^-iE;;t;;\//£j-^ forz = r + l,---,r + s ' ^'^ 

(1 2) 

In deriving the above expressions for A^j, we have used the relationships between Qg ' 
and Qo"* similar to those given in (2.13). Keep in mind that q} in (3-1) as realized by 
(2.15) (and (2.16)) form the A/" + 1 dimensional representation of the polynomial algebra 
(2.17). This representation is also realized by the differential operators (2.25) acting on 
the -I- 1 dimensional space of polynomials with basis |l, 2;, 2;^, We can thus 
equivalently represent (3.1) (i.e. (1.1)) as the single-variable differential operator of order 
M = max{E[=iA:^,E[;;+iA;„2}, 

r+s r+s 

+9^ ^ n , n % (-^ - - 1 + 2" + ^? - ''"'i}'^' ) 
+j.-'nn%f4+,,+.«-(i^ii^) (3.3) 



with 



j=ii=i \ f^j 



h{zj- + qr + sf^)-l fori = l. 



[ ki{2K - - Qr - t + si ') - for i = r + 1, • • • , r + s 
We will now solve for the Hamiltonian equation 

Hijjiz) ^ EiP(z) (3.5) 

by using the Functional Bethe Ansatz method, where ip{z) is the eigenfunction and E is 
the corresponding eigenvalue. It is easy to verify 

Hz"^ = z^+'g n f] kj (-n-qr-t + 2n + sf - ^l'^' + ^ ] 

j=r+li=l \ % / 

-I- lower order terms, n e Z+. (3.6) 
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This means that the differential operator (3.3) is not exactly solvable. However, it is quasi 
exactly solvable, since it has an invariant polynomial subspace of degree jV + 1: 



HVCV, V = span{l,z,...,z^}, dimV = J\f + 1. (3.7) 

This is easily seen from the fact that when n — jV, the first term on the r.h.s. of (3.6) 
becomes z^^^ g YTjt^+i IliLi kj (^Qr+s + sf^ — ■^^n^?^^ which vanishes identically for all 

the allowed g^+s values (noting sf.^^ = 0). 

As (3.3) is a quasi exactly solvable differential operator preserving V, up to an overall 
factor, its eigenfunctions have the form, 

AT 

ij{z) = l[{z-a,), (3.8) 

1=1 

where {ai\i ~ 1,2, •• • ,J\f} are roots of the polynomial which will be specified later by the 
associated Bethe ansatz equations (3.11) below. We can rewrite the Hamiltonian (3.3) as 

/ dV 

H = Y.P^i^)[-^^) +Poiz), (3.9) 

where Pq{z) and Pi{z) are polynomials in z resulted from the expansion of the products 
in (3.3). 

Dividing the Hamiltonian equation Hip — Eip over by ip gives us 
Hil) MM 1 

E-^ = ZP>(-)« E ) (3.10) 

The l.h.s. of (3.10) is a constant, while the r.h.s is a meromorphic function in z with at 
most simple poles. For them to be equal, we need to eliminate all singularities on the 
r.h.s of (3.10). We may achieve this by demanding that the residues of the simple poles, 
z = ai.i = 1,2, ...,J\f should all vanish. This leads to the Bethe ansatz equations for the 
roots {cti} : 

E E 7 7 _ r + Pi(ap) = 0, p^l,2,...,M. 

(3.11) 

The wavefunction ip{z) (3.8) becomes the eigenfunction of H (3.3) in the space V provided 
that the roots {ai} of the polynomial ■0(2;) (3.8) are the solutions of (3.11). 

Let us remark that the Bethe ansatz equation (3.11) is the necessary and sufficient 
condition for the r.h.s. of (3.10) to be independent of z. This is because when (3.11) is 
satisfied the r.h.s. of (3.10) is analytic everywhere in the complex plane (including points 
at infinity) and thus must be a constant by the Liouville theorem. 
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To get the corresponding eigenvalue E, we consider the leading order expansion of 



It is easy to show that V±fiip{z) have the expansions, 



-n 



j^ LLp=r+l '"p 



X 



n n 

j=r+l 1=1 \ 



V-tp — z 



1 + s 



(2) 



N 
1=1 



r+s 

S 



-nK+-^-+'-f'- ''"'.'fi"" i + 



j=r+l 



(3.12) 



Substituting these expressions into the Hamiltonian equation (3.5) and equating the 
terms, we arrive at 

r ^ 1 \ ^ '^"'"^ / 1 \ ^ 

E = Y.wu(ki{N + qr + sf'^)--) + E WijA;i(g,+, + sf^) - - 

i=l ^ '^i^ i=r+l ^ '^i^ 



i<j=2 



+ E ( ^^l-A/" + ^r- + sf^) -i:)(kjiAf + qr + Sf^) - 



+ E ^ij {HQt+s + sf^) - ( fci(gr+. + ^) - ) 

i<j=r+2 ^ /^i/ V '^3 J 

+ E «^ J ^^(•A^ + ?r + S?^) - TT + E ^4 ^^(?-+- + «?^) - 
i=l ^ i=r+l ^ 



j=r+l 1=1 



^1 



E«i' 

i=l 



(3.13) 



where {ctj} satisfy the Bethe ansatz equations (3.11). This gives the eigenvalue of the 
2-mode boson Hamiltonian (1.1) with the corresponding eigenfunction ^l){z) (3.8). 



4 Explicit examples corresponding to BECs 

In this section we give explicit results on the Bethe ansatz equations and energy eigen- 
values of the Hamiltonian (1.1) for the special cases of r + s < 4. These models arise 
in the description of Josephson tunneling effects and atom-molecule conversion processes 
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in the context of BECs. Our approach provides a unified treatment of their exact solutions. 



A. r = 2, s = 1, ki — k2 — ks — 1 
The Hamiltonian is 

3 3 

H '^'■i.^'-i. + WijNiNj + g (a| + aia24) • (4. 1) 

i i<j 

This is the so-called hetero- atom-molecule BEG model and has been solved in [23, 24] via 
the Algebraic Bethe Ansatz method. Here for completeness we present the exact solution 
derived from our approach without providing any details. 
The Bethe ansatz equations are given by 

+ + ,.i,2,.,A. (4.2) 

i^pOii-ap gap + Aua^ 

with A/" = 2fi; — 2 — Y and the energy eigenvalues are 

E = Wll{N+llf + w22^f^ + Wl{^^+h) + w2^^+Wl2^^{N+h)-gYa^. 



Here 



(4.3) 



^11 = W22 + Wll + W33 W12 - Wis - ^23, 

Bii = W1-W3 + W2 + W22 + wn (2/1 + 1) + W33 (5 + /i - Ak) . (4.4) 



B. r = 2,s = l, ki = k2 = l,k3 = 2 
The Hamiltonian is 

3 3 

H^J2 '^i^i + E '^ij^i^j + 9 {aiaUl + aia2af) . (4.5) 

i i<j 

This yields a model of three-mode atomic-molecular BECs which has not been exactly 
solved previously. Specializing the general results in the preceding section to this case, 
we have ^2 = 1, Q's = 4 or |, t = ^ = y and J\f = 2k — 1 — — The differential 
operator representation of the Hamiltonian (4.5) is thus 

where 

P2{z) = 4gz^ + A2iz^ + gz, 

Pi{z) = S2iz2 + L»2iz + ^(Zi + l), 

P^{z) = F21Z (4.7) 
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with 



^21 = Wn + W22 - 2^23 + 4^33 - 2^13 + Wi2, 

B21 = %(4 + /i-4k), 

D21 = W2 + W1+ W22 + wu {h + 1) - 2^3 + W33 (14 + 4/1 - 16k) 
+wu {2h + 1) + Wis - ^ - 3^i) + W23 (ak-^- , 

F21 = g(^{4K-2-h){4K-4-h) + ^^ 



G21 = wih + W3 ^4k ~ 2 ~ ^0 ^^^^^ (^'^ 2 

+Wn/2 + W33 (^(4/^ - /i - - j . (4.8) 
The Bethe ansatz equations are given by 

E^^ ^f tfrtf P^l,2,-..,A. (4.9) 

^ a, - ftp 45^0;^ + ^21 ftp + 

and the energy eigenvalues are given by 

/ 1\^ 

+2 [q-i - ^ [wi3 [M + /i) + W2-iM + W3] + wi (A^ + /i) + W2H 
-4^(?3 + J) (93 + ^)E«- (4-10) 



where A/" = 2k - 1 - q'3 - ^. 
C. r = s = 2, /ci = A;2 = ^3 = ^4 = 1 
The Hamiltonian is 

4 4 

if = E WiNi + E WijNiNj + g (a\ala3a4 + aia244) . (4-11) 

i i<j 

This gives a model of four-mode atom-molecule BECs. This model has not been exactly 
solved previously. Applying the results in the preceding section gives q2 = ^ — Qa, 
t = lis^i^ + sf^) = |(/i + k) and A/" = 2/s: - 2 - i(/i + /3). The differential operator 
representation of the Hamiltonian (4.11) thus reads 

H = P2{z)^ + Pi{z)^ + Po{z), (4.12) 

where 

P2{Z) = g{z^ + Z) + A22Z\ 

Pi{z) = B22Z'' + D22Z + g{li + l), 

Po(z) = gX(J\f + h)z + G22 (4.13) 
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with 

^22 = Wii + W33 + W22 - W24 - Wi3 + W12 - W23 + W44 - Wu + W34, 

B22 = g{h + 5-AK), 

D22 = wi + tii2 - - + W22 + W33 (1 - 2/3 - 2A/') + wii (2/1 + 1) 
+W12 (1 + h) + W23 (A/" + h - 1) + (A/" + h - /i - 1) 
+^14 (TV - - 1) + W44 (1 - 2AA) + 'u;34 (1 - ^3 - 27V) + W2aM, 

G22 = Wih + (TV + Zg) + Wj^ + (TV + ^3)^ + W^^f W + ^3) 

+^11/1^ + WuhN + ^447^^ + -^^13^1 (TV + h) , (4.14) 

where (and below) j'V = 2k — 2 — |(/i + ^3). 
The Bethe ansatz equations read 

^ 2 ^22Q^^ + /^22«p + ^7(/l + l) ^_ 

= ~2 ' P=1>2,---,A/ (4.15 

^ tti - ftp 5((a^ + ftp) + A22a^ 

and the energy eigenvalues are 

E = 'W;ii(7V + Zi)V'U;227V^ + 'U;33^3+'«^127V(7V+^l)+'»^13^3(7V+Zi) 

+ {w2zh + W2)N + wi {Af + h) +W3I3- g{h + 1) I] ttj- (4.16) 



i=l 



5 Discussions 

We have introduced the higher order polynomial algebra (2.17) and identified it as the 
dynamical symmetry algebra of the multi-boson system (1.1). We have constructed its 
unitary irreducible representation (2.23) and the corresponding single- variable diff'erential 
realization (2.25). We have then used the differential realization to rewrite the Hamilto- 
nian (1.1) as a QES differential operator (3.3) acting on the finite dimensional monomial 
space, thus providing an algebraization of the higher order Hamiltonian differential equa- 
tion (3.5). Exact eigenfunctions and eigenvalues of the Hamiltonian (1.1) have been found 
by employing the Functional Bethe Ansatz technique. As examples, we have provided ex- 
phcit expressions for the BEC models which correspond to the r -|- s < 4 cases of (1.1). 

It would be interesting to adapt the techniques of this paper to obtain exact solutions 
of the generalized Lipkin-Meshkov- Click and Tavis-Cummings type model defined by 

H = ^w,iV, + f]Jo^ + ^(4ati---a^ + J^a^-.-aJ'^'^), (5.1) 

i=l 

where J±,o are the generators of the su{2) algebra. Also we plan to explore the q-boson 
counterparts of the multi-boson systems and their underlying q-deformed polynomial 
algebraic structures. Lastly, it is of interest to further explore the connection between the 
higher order ODEs and integrable models (i.e. the ODE/IM correspondence) along the 
fine of [28]. Results in those directions will be presented elsewhere. 
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